Complete Determination of the Spectrum of a Transfer Operator associated with 

Intermittency 



O 

o 

(N 

< 

(N 



U 



> 
o 

00 

O 

o 

"a 

> 
X 



Thomas PrellbergQ 
Institut fiir Theoretische Physik, Technische Universitdt Clausthal, 
Arnold- Sommerf eld Strafle 6, D-38678 Clausthal-Zellerfeld, Germany 
(Dated: February 8, 2008) 

It is well established that the physical phenomenon of intermittency can be investigated via 
the spectral analysis of a transfer operator associated with the dynamics of an interval map with 
indifferent fixed point. We present here for the first time a complete spectral analysis for an example 
of such an intermittent map, the Farey map. We give a simple proof that the transfer operator is self- 
adjoint on a suitably defined Hilbert space and show that its spectrum decomposes into a continuous 
part (the interval [0, 1]) and isolated eigenvalues of finite multiplicity. Using a suitable first-return 
map, we present a highly efficient numerical method for the determination of all the eigenvalues, 
including the ones embedded in the continuous spectrum. 

PACS numbers: 02.30.Sa, 02.30.Tb, 02.70.Hm, 05.10.-a, 05. 45. Ac 



INTRODUCTION 

Intermittency, one of the main routes from order to 
chaos Q, is characterized by the loss of stability of a 
fixed point of the dynamics. The dynamics directly at the 
transition is determined by a marginally unstable fixed 
point near which trajectories are slowed down severely, 
leading to the for intermittency characteristic interplay 
of chaotic and regular dynamics. Such a behavior can 
easily be modeled by a map / of the unit interval [0, 1] 
which is uniformly expanding everywhere except near an 
indifferent fixed point at zero, where f{x) ^ x + cx^~^^ 
as a; ^ with exponent 1 + ?' > 1 . A typical example is 
the Manneville map f{x) = x + x i+'-modl II . 

There have been many theoretical approaches to the 
description of intermittency such as renormalization 
group analysis Q . An approach suited to rigorous treat- 
ment is given by the thermodynamic formalism and 
leads to the spectral analysis of transfer operators. In 
contrast with uniformly expanding maps, however, the 
indifferent fixed point induces non-Gibbsian equilibrium 
states or, more precisely, weakly Gibbsian states [||. 
Therefore, deeper understanding based on rigorous anal- 
ysis has been hard to come by. Some progress was made 
by studying a suitably defined piecewise linear interval 
map Q, albeit at the cost of simplifying the dynamics 
by severely reducing correlations. In 0, || we argued the 
disappearance of a spectral gap for the Perron-Frobenius 
operator, implying loss of an exponential decay of corre- 
lations for the dynamics. The decay of correlations has 
later been shown to follow a power law; in the analytic 
case (r = 1) a numerical estimate of has been 

complemented by a rigorous upper bound of t^^logt, 
obtained by random perturbation techniques Only 
very recently has the Perron-Frobenius operator for a 
particular intermittent map, the Farey map, been shown 
to be self-adjoint on an appropriate Hilbert space Q 
with continuous spectrum on the interval [0,1]. There 



is also numerical work available describing just how the 
continuous spectrum emerges when approaching the in- 
termittency transition p^ . 

Much less is known about the spectrum of the Ruelle- 
Perron-Frobenius (RPF) operator, which generalizes the 
Perron-Frobenius operator within the framework of the 
thermodynamic formalism. In ^ we showed that for a 
general class of intermittent maps this operator is quasi- 
compact with essential spectral radius equal to one, and 
that the leading eigenvalue undergoes a phase transition 
characteristic for intermittent dynamics. Only recently 
was it shown that for a class of piecewise analytic maps 
the continuous spectrum is restricted to the interval [0, 1] 
on an appropriate function space p^ . The present work 
gives for the first time a complete spectral analysis of the 
RPF operator for an intermittent map, the Farey map. 

A promising conceptual approach to the study of in- 
termittency is the study of a first-return map (or induced 
map) with respect to a domain of phase space away from 
the intermittent region. In a pioneering study ^ we 
introduced a suitably modified transfer operator for the 
induced map and related its spectral properties to the 
ones of the transfer operator of the intermittent system. 
This approach has since been extended to the study of 
other quantities such as regularized Frcdholm determi- 
nants and dynamical zeta functions 14, |l^, ^ and 
multi-dimensional systems One essential idea em- 

ployed is that one can understand a general transfer op- 
erator V of an intermittent map by first considering its 
"intermittent part" Vq and then viewing the "chaotic re- 
mainder" 7^1 = 1 — T' as a perturbation. In this way, 
many of the results presented here can in principle be 
extended to more general intermittent maps, although 
we shall focus our attention on the Farey map of the in- 
terval [0, 1] onto itself, which is defined as 



/(^) 



foix) 



x/{l-x) 
(1 - x)/x 



if < a; < 1/2, 
if 1/2 <x< 1. 



(1) 



2 



We denote the inverses by Fo{x) ^ fo ^{x) = x/{l + x) 
and Fiix) = fr'^ix) = 1/(1 + a;). 

The advantage of this map is that each branch can 
be analyticaUy extended to the whole Riemann sphere 
and that higher iterates of the left branch can be given 
exactly, fo^{x) = x/{l — nx), a fact that is intimately 
connected to the global conjugacy of fo{x) to the shift 
X ^ X — 1 p3| . The thermodynamic formalism suggests 
the study of the RPF operator P associated with a map 
/, which is given by Vip{x) = J2f{y)=x \f'iy)\~'^'Pifi^)) 
with /9 G K. For the Farey map the operator consists 
of two terms which can be readily identified with the 
"intermittent" and "chaotic" parts. We thus write V = 
Vo+Vi where V^lp{x) = \F,' {x)fip{F,{x)) and ^ G R, 
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Despite the apparent simplicity of this operator, its 
properties have not been well understood until recently. 
In this paper we first present an abstract analysis de- 
scribing its spectrum and then proceed to determine all 
of its eigenvalues with a highly efficient numerical method 
based on the first-return map. 



SPECTRAL PROPERTIES OF THE TRANSFER 
OPERATOR 

As indicated above, one can understand the spectral 
properties of V by first studying the operator Vq and 
then viewing Vi as a perturbation. 

Under conjugacy with Cip{x) — ip{fi{x)), the operator 
Vo transforms to the shift operator Sip{x) = ip{l + x). 
This operator is now independent of /3, and on a suitably 
defined function space has continuous spectrum a{S) = 
(Tc{S) = [0,1]. This can be understood by considering 
functions ip which are obtained by a generalized Laplace 
transform (p{x) = C'ijj{x) = exp(— sa;)V'(s)(i/i(s) of 
square integrable functions ip G L'^{R+,fi). (This au- 
tomatically ensures analyticity of ip{x) in a suitable do- 
main.) The action of the operator S is then conjugate 
to multiplication by exp(— s) on L'^ {M.+ , jj.) . The spec- 
trum of this multiplication operator is continuous and 
given by the closure of the range of the multiplying func- 
tion. This proves that Vq is a bounded self-adjoint op- 
erator on the Hilbert space C CL'^ , ji) with spectrum 

To study the spectrum of 7^, we consider the identity 



l-zV = il~zro)il-M,) 
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with Aiz = (1 — zVo)^^ zVi being an operator- valued 
analytic function for z G C — [1,cxd) We point out 
that (|^) remains valid when A^^ is analytically contin- 
ued across the cut [1, oo). Expanding formally in powers 



of z, we find that Mz = J2n=i z'^Va^-^Vi is a transfer 
operator associated with the induced map g on [1/2, 1] 
with the branches g„ = /o"^^/i for ri G N, where the re- 
turn time n has been encoded via a multiplicative weight 
factor z". 

To determine the spectral properties of the operator 
M.Z we note that the induced map g on [1/2,1] is ex- 
panding. It follows that for \z\ < 1 the operator Aiz 
acting on a Frechet space of functions analytic in a do- 
main containing the interval [1/2, 1] is nuclear of order 
zero. Therefore, the analytic continuation of A4z is also 
nuclear of order zero (which can be seen via the analytic 
continuation of the associated Fredholm determinants, 
seell). 

The identity (^) immediately implies that X = z ^ is 
an eigenvalue of V if and only if 1 is an eigenvalue of ^ . 
The respective eigenspaces are identical, so that the ge- 
ometric multiplicities of the respective eigenvalues z~^ 
and 1 are the same. Using the analytic continuation of 
Aiz across the cut [1,cxd), general analyticity arguments 
and the nuclearity of imply that the non-zero point 
spectrum of V consists of isolated eigenvalues with finite 
multiplicity, with and 1 as the only possible accumula- 
tion points. Moreover, bounds on A4z for jzj < 1 imply 
that there are only finitely many eigenvalues with mod- 
ulus greater than one. Along the cut z G (l,oo) we see 
that if the analytic continuation of does not have 
an eigenvalue 1 then 1 — Mz is invertible, implying that 
cTciV) = tTc(7^o) = [Ojl]; albeit with the possibifity of 
embedded eigenvalues. Similar arguments can be found 
in a rather general setting in [^ . We further know that 
zero is an eigenvalue of infinite multiplicity, as Vtp — 
for any (p satisfying ip{x) — —(p{l — x). 

We now show that V is in fact self-adjoint on a suit- 
ably defined Hilbert space for which we can give a rep- 
resentation of V which is explicitly symmetric. For 
this, we consider the transfer operator Q for a gen- 
eral Mobius transformation h(x) — (ax + h)/{cx + d) 
with a,b,c,d G M and a = ad — be = ±1. We have 
Q(p{x) — {a — cx)^^^ Lp{^^^) and under a generalized 
Laplace transform Q is conjugated to an integral oper- 
ator /CV'(s) = K{s,t)^P{t)diJL{t) onL2(R+,^). Using 
the Schlafii integral representation for Bessel functions 
||l8| , one can calculate for c 7^ directly 
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with Z^{u) — Ii,{u) for a = 1 and Z^{u) ~ J^{u) for a = 
— 1. Demanding that the kernel be symmetric determines 
the appropriate measure. In this way, we obtain for Vi 
the measure /i(s) — e~'^s'^^~^ and the kernel K{s,t) — 
{st)2~l^ J2i3_i[2y/si). Combining this with our discussion 
of the operator Vq^ we get a conjugacy of the transfer 
operator V to the self-adjoint operator 

/>oo 

5V^(s) = e"^V(s) + / K{s,t)ij{t)dfi{t) (5) 
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Expressed differently, V is self-adjoint on the Hilbert 
space 

. 

with /i(s) = e ^s^'^ ^ and appropriate induced inner 
product. M.Z also leaves this Hilbert space invariant, 
and via a similar argument one can show self-adjointness 
for real z <\. 

The case (3 = 1 has already been treated in al- 
though with a slightly different choice of measure, mo- 
tivated by the Hilbert space approach for the continued 
fraction transform There, it was shown that 

for /2(s) = s/(exps — I) the only non-zero eigenvalue 
of 7^ is I with eigenfunction ip{x) = 1/x. The measure 
fi{s) = se'" considered here has the advantage that the 
self-adjointness is explicitly evident from (||). However, 
f{x) = 1/x is not an element of the Hilbert space con- 
sidered here, as it corresponds to 't(i{s) = 1/s which has 
infinite norm in L^(IR_|_,^). We find this quite natural, 
as the usual interpretation of this eigenfunction is as an 
absolutely continuous invariant density with respect to 
the Lebesgue measure on the interval [0, 1], which in this 
case is non-normalizable. 

NUMERICAL ANALYSIS OF THE SPECTRUM 
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FIG. 1: The seven leading eigenvalues of Mz for (5 = 0.5. 
The analytic continuation of the eigenvalue branches beyond 
z = 1 is complex valued and has a jump across the cut [1, co). 
Only the largest eigenvalue branch intersects A = 1. 
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As shown above, Mz is an operator- valued analytic 
function with nuclear spectrum. By standard analytic 
perturbation theory ||2^, the eigenvalues of M.^ are 
(branches of) analytic functions in z with only alge- 
braic singularities. This therefore provides a possibility 
to compute the eigenvalues of V numerically; choosing 
an eigenvalue branch A„(z), one simply needs to solve 
A„(z) = 1 to get an eigenvalue A = 1/z oiV. The essen- 
tial advantage over working directly with V is that we 
have removed the continuous spectrum which presents 
enormous difficulties for a direct numerical analysis. 

As the operator Mz acts on a Hilbert space of analytic 
functions, it is reasonable to consider the action of A^z on 
coefficients of power series. From the explicit expansion 
MzVix) = Er=i^"(l + nx)-''^^(l - x/(l + nx)) we 
obtain by expanding Lp{x) in a power series around x = 
1 matrix elements A^"'™ in terms of the polylogarithm 
Li.(^) =E^=i|^ as 




We then approximate Mz with truncated operators 
A^i^'' acting on a subspace of polynomials of at most 
degree N . It turns that this approximation works very 
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FIG. 2: The seven leading eigenvalues oi Mz for /3 = — 3. 
Here, the eigenvalue branches can be continued beyond 2 = 1 
and intersect A = 1 at z = -2.971, z = -0.168, z = 0.038, 
and z — 13.101, the last value being outside the range of this 
plot. 



well and that one obtains the values of the leading eigen- 
values to a high accuracy. In this setting, this idea goes 
back to |l^. In where a related operator was stud- 
ied, the leading eigenvalues have been obtained in this 
way with accuracy of 10^^^ p^ . 

We begin the description of the results of our analysis 
by briefly considering the special cases (3 = —N/2 with 
N G Nq. For these values of /3 we find that polynomials of 
at most degree N give an {N + l)-dimensional invariant 
subspace for Mz and P. The corresponding truncated 
matrix A^l^^ has entries which are rational functions in z 
from which one can easily calculate N +1 eigenvalues ex- 
actly. Numerically we observe that eigenf unctions oi Mz 
which are not in this invariant subspace have eigenval- 
ues which are strictly smaller in modulus. It is especially 
noteworthy that these leading eigenvalues can be analyt- 
ically continued across z = 1, which allows for eigenval- 
ues of the Farey operator V embedded in the continuous 
spectrum. However, a more detailed numerical analysis 
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terplay of sub-dominant eigenvalues and the continuous 
spectrum. Both the rigorous arguments and numerical 
methods are generalizable to more complicated systems 
with interniittency. 



FIG. 3: The spectrum of P as a function of 13, showing 
two non-zero eigenvalue branches and isolated eigenvalues at 
/3 = — 1 and (3 = —2 which are embedded in the continuous 
spectrum [0, 1]. 



shows that this analytic continuation is possible only for 
/3 = —N/2, and a small deviation from these values leads 
to eigenvalue branches A„(z) whose analytic extension 
generically shows a non- vanishing jump in the imaginary 
part along the cut [l,oo) (26|. Therefore \n{z) = 1 for 
real z > 1 can only be satisfied when (3 — —N/2, and no 
eigenvalues of V embedded in the continuous spectrum 
exist for other values of (3. Thus, in what follows we shall 
mainly be concerned with real z < 1. 

For /3 > 1 we find that all the eigenvalues oi M.z are 
strictly less than one in modulus, so that V has contin- 
uous spectrum [0, 1] and no non-zero eigenvalues at all. 
As (3 decreases the eigenvalues oi increase in value. 
For /3 > 1, we find that the leading eigenvalue branch 
oi M.Z intersects A = 1, leading to the emergence of a 
simple leading eigenvalue of V for (3 < 1 The z- 

dependence of the spectrum of A^z is shown in Figure 
|l| for /3 = 1/2. Only the largest eigenvalue intersects 
A = 1, implying that V has only one non-zero eigen- 
value. At (3 — —3/2 a second eigenvalue branch begins 
to intersect A = 1 at large negative z, implying that V 
has two non-zero eigenvalues. The second eigenvalue of 
V is negative and becomes in modulus equal to the essen- 
tial spectral radius aX (3 — —2, and therefore determines 
for (3 < —2 the spectral gap which controls the decay of 
correlations. Figure ^ illustrates the leading spectrum of 
for (3 — —i. At this special value of /3 we can extend 
the eigenvalue branches beyond z = \ and find a total of 
four eigenvalue branches which cross A = 1. One of these 
crossings is at z = 13.101, corresponding to an eigenvalue 
of V embedded in the continuous spectrum. Proceeding 
in this fashion, we can numerically determine the com- 
plete spectrum of V for arbitrary real [3. Figure ^ shows 
the spectrum obtained in this way for —2.5 < /3 < 1.5. 

In summary, we have abstractly characterized the spec- 
trum of the transfer operator for the Farcy map and pre- 
sented a highly efficient method for the explicit compu- 
tation of its eigenvalues. Of special interest is the in- 
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For general parabolic fixed points such a conjugacy can 
always be found locally, a fact that is used in the analysis 
of pi- 

In 111] a related identity has been used, leading to the 
study of A4i = Pi(l - zPo)"^ Obviously, Mz and M'^ 
are related by conjugacy. 

As in in [^^, we find that we can improve convergence 
by expanding around a different point. However, this 
leads to the introduction of spurious eigenvalues which 
make the analysis more difficult. By choosing an expan- 
sion around a; = 1 we avoid spurious eigenvalues and find 
that, overall, the spectrum is better approximated. 
One can understand this by considering the analytic ex- 
tension of \As{z), which across the cut [1,00) jumps at 
2; > 1 by an amount 27ri log""^ (2:)/r(s). 
In 1^, g] it was proved that this leading eigenvalue de- 
creases to 1 like —(1 — /3)/log(l — /3) as /? approaches 1 
from below. 



